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Abstract 

We investigate multiple Charlier polynomials and in particular we will use the 
(nearest neighbor) recurrence relation to find the asymptotic behavior of the ratio of 
two multiple Charlier polynomials. This result is then used to obtain the asymptotic 
distribution of the zeros, which is uniform on an interval. We also deal with the 
case where one of the parameters of the various Poisson distributions depend on the 
degree of the polynomial, in which case we obtain another asymptotic distribution 
of the zeros. 

1 Introduction 

Charlier polynomials {C^\n = 0,1,2,...} are orthogonal polynomials for the Poisson 
distribution, i.e., 

E^W^W^f = 0, n^m, 

k=0 

where a > 0. These polynomials are orthogonal on the positive integers and as a result 
their zeros are separated by the integers: between two consecutive integers there can 
be at most one zero of Cn . Charlier polynomials have various applications, e.g., in 
queueing theory [TT] and recently [23] , in the analysis of the lengths of weakly increasing 
subsequences of random words |10j . and in the totally asymmetric simple exclusion process 
(TASEP) [2j. Their asymptotic behavior has been studied by Maejima and Van Assche 
[14] . Kuijlaars and Van Assche [12] . Rui and Wong [20] . Goh [7J, Dunster [5] and most 
recently by Ou and Wong [IT] using the Riemann-Hilbert method. 

We will investigate multiple Charlier polynomials, which are polynomials of one vari- 
able with orthogonality properties with respect to more than one Poisson distribution. 
Take r Poisson distributions with parameters a%, . . . , a r > and such that cij ^ a,j when- 
ever i j. Let n = (n\, n 2 , . . . , n r ) be a multi-index of size \n\ = n\ + n 2 + • • • + n r , then 
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the multiple Charlier polynomial is the monic polynomial of degree |n| for which ([TJ 
p. 29-32], p. 632], [2]]) 



J2 c ^ k ~ir = i 1=0,1,-^-1, j = i 

k=0 



r. 



For r = 1 we retrieve the Charlier polynomials. The multiple Charlier polynomials can 
be obtained using the Rodrigues formula [TJ |9j [21] 



Cn(x) = (-1)1* ( Ha^j T(x + 1) ^n^"^ * ) ffr + i) 



where V is the backward difference operator, given by Vf(x) = f(x) — f{x — 1). An 
explicit formula for the multiple Charlier polynomials is 



C rt(x) = " ' ^(- n l)ki ■ ■ ■ (-nr)k r (-x) kl +k 2 +-+k r 

k 1= k r =0 

(-a 1 ) ni - kl (-a 2 ) n2 - k2 ■ ■ ■ (-a r ) nr 



k l \k 2 \---k r \ 



;i.2) 



Multiple Charlier polynomials satisfy a number of (higher order) difference equations 
(Lee [13] and Van Assche [21]). They appear in remainder Pade approximation for the 
exponential function [Tj5], as common eigenstates of a set of r non-Hermitian oscillator 
Hamiltonians [T5], and we believe that they are related to the orthogonal functions ap- 
pearing in two speed TASEP (totally asymmetric simple exclusion process) [3]. 

In this paper we first obtain in Section 2 some properties of the multiple Charlier poly- 
nomials, such as the generating function and the nearest neighbor recurrence relations. 
The zeros of multiple Charlier polynomials are real, positive and separated by the positive 
integers, as is the case for the usual Charlier polynomials: between two positive integers, 
there can be at most one zero of a multiple Charlier polynomial (see, e.g., [T5| Theorem 
3.4]). The largest zero of is therefore > \n\ — 1. In order to prevent the zeros to go to 
infinity, we will use a scaling and consider the scaled polynomials Pu,n(x) = Cg(Nx)/N^. 
One of the main results in this paper is in Section 3 where we obtain the asymptotic be- 
havior of the ratio of two scaled neighboring multiple Charlier polynomials. We use that 
result in Section 4 to obtain the asymptotic zero distribution of the scaled multiple Char- 
lier polynomials. Another important result is in Secton 5 where we give the asymptotic 
behavior (ratio asymptotics and zero distribution) when one of the parameters depends 
on the scaling N. This gives a different asymptotic zero distribution which is somewhat 
more interesting. 

2 Some properties of multiple Charlier polynomials 

2.1 Generating function 

Charlier polynomials have the generating function [U Ch. VI, Eq. (1.1)] 

00 j-n 

^ C '" a)(X) d =(1+tre " at ' |t|<L (2 ' 1} 

n=0 
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For multiple Charlier polynomials one has a multivariate generating function (with r 
variables) . 

Theorem 2.1. Multiple Charlier polynomials have the following (multivariate) generating 
function 



,ni.n 2 j.n r 
■ x \ 1 2 " " " l r 

' ni!n 2 ! • • -n r \ 

rii=0ri2=0 n r =0 



EE-la 



= (1 + *! + t 2 H 1- t r ) x exp(-ai*i - a 2 t 2 Mr)- ( 2 -2) 

Proof. We can use induction on r. For r = 1 we have the familiar generating function for 
Charlier polynomials ( 12. II) . 

Suppose the result is true for r — 1, then observe that 111. 2ft implies 



-a,. 



k r =0 



Hence the multivariate generating function is 
(1 + t x H h t r -x) x exp(-ai*i a r -xU-\) 



£ c 1 + *i + • • • + ^r^u-^k^^ 

n r =0 fc r =0 



Changing the order of summation gives 
(1 + t\ H h tr-i) x exp(— aiti - • • ■ - a r _ii r _i 



\ 00 j-n 



and by putting i = n r — k r 



(1 + tx H h t r -i) a exp(— aiti - ■ ■ ■ - a^-iV-i) 

Ed +*. + •■■ +t^.)^(. I )f; 

fc r =0 ^ r ' £=0 

Now use 



and 



S (1 + * + " ■ + ^ £) - + i +tl+ ^ + ... +<r J ; 



Av=0 

to obtain the desired result. □ 
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The region of convergence of this generating function is a log-convex set in C r , which is 
the case of all power series in several variables, and the series certainly converges whenever 
\tj\ < 1/r for every j 6 {1, 2, . . . , r}, or when \tj\ < Cj for 1 < j ' < r, where < Cj < 1 
and Ej=i c j — 1- As a corollary, one can obtain an integral representation of the multiple 
Charlier polynomial, by integrating r times over a closed curve around 0: 

Ch(x) 1 £ I (1 + Z\ H V z r ) x exp(-aizi - ■ • • - a r z r ) 

dz\ ' ' ' dz r . 



m\---n r \ (2ni) r J J z^ 1+1 ■ ■ ■ z n - +1 

2.2 Recurrence relations 

For multiple orthogonal polynomials there is always a nearest neighbor recurrence relation 
of the form 

r 

xPh(x) = Pn +e - k {x) + bfi jk Pfi(x) + a H jP H _ gj (x) (2.3) 

where k — l,...,r Thm. 23.1.11], [22], and e k = (0, . . . , 0, 1, 0, . . . , 0) is the kth unit 
vector in N r . The recurrence relation for multiple Charlier polynomials was given in [HI 
p. 632] without proof. Here we will work out the details of the proof. 

Theorem 2.2. The nearest neighbor recurrence relation for multiple Charlier polynomials 
is 

r 

xC H (x) = C H+Sk {x) + (a k + \n\)C H {x) + ^ n j a j Cf t s.{x). (2.4) 

Proof. From flT2} and (-x) n = (-l) n x n + (-1)™" 1 (^x"" 1 + • • • we find that 

C n (x) = xW + SnxW- 1 + ■ ■ • , 

where 5^ can be found by taking /c 2 , . . . , k r ) = (n\, n 2 , . . . , n r ), which gives the con- 
tribution — ('2') t° fifi, and for each j with 1 < j < r we get for (fc 1; k 2 , . . . , k r ) = 
(rii, n 2 , . . . , Uj — 1, . . . , n r ) the contribution —ajUj, so that 



r 



If we compare the coefficient of x' n ' in (I2.3p . then bft,k = <^n ~ <5n+e fc , which for the multiple 
Charlier polynomials gives b^ — \n\ + For the recurrence coefficients a^j we can use 
Eq. (23.1.23)] 

Er=ofc WJ ^(fc)a 3 Vfc! 

The sums can be computed using the Rodrigues formula (11.11) : the difference operators 
a~ x V ni a x (i — 1, 2, . . . , r) are commuting, so we can first apply a~ x V nj a x to find 

00 a k r oo / r \ 

E*" ,c, *(*)^ = (- 1 ) w II o i < E A:BJ ( VBJa i) II a * rfevn<a *hr 

fc=0 ' i=l fc=0 \i=l,i^j / 
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Now use summation by parts times to find 



oo k r oo / r \ 

T.wcnW^ = (-i) |ni n a r(-i) nj E( Anj ^>? n °r fc v n x fe U 

fe=0 ' ' i=l k=0 \i=l,i/j / 

r oo / r \ 1 

= (-i) w nc(-wE fl * n a ^ vn ^hi- 

i=l k=0 \i=l i^j / 



If we change n,- to n,- — 1 then this gives 



fc=0 



k 

k\ 



i-ipmci^i-ir 1 ^-!)!!;^ n a ^ a ^\lv 

\i=l J k=0 \i=l,ijLj / 



Dividing both expressions then gives 



□ 

The recurrence coefficients are quite simple in this case, and in particular a^j = njcij > 
whenever rij G N. This implies that the zeros of Ca and its nearest neighbors C r t+e fe 
interlace for every k G {1, 2, ... , r}, see [8]. This will be useful in the next section. 

3 Ratio asymptotics 

There are various levels of asymptotic behavior to consider. In this paper we limit the 
analysis to ratio asymptotic behavior, i.e., the asymptotic behavior of the ratio of two 
neighboring polynomials. In order to prevent the zeros from going to infinity, we use a 
scaling and we will investigate the ratio Cft + g k (Nx) /Ca(Nx) for x G C \ [0, oo), where N 
is of the order \n\, i.e., limAr^oo \n\/N = t > 0. 

Theorem 3.1. Suppose rij = \_qjn\, with < qj < 1 and Yl^iQj = 1? so ^at \n\/n — > 1 
as n — >■ oo. Let a« > for 1 < i < r and ^ aj whenever i ^ j . Then for t > and for 
every k G {1, 2, . . . , r} one has 

Jim C n + e k {Nx) =x _ t 
n^oo,n/N^t NCh(Nx) 

uniformly for x G K , where K is a compact set in C \ [0, oo). 

Proof. We will use the notation Ph,n(%) — Cft(Nx)/N^ for the monic and rescaled mul- 
tiple Charlier polynomials. The zeros of Cfi~s- and are real, positive and interlace 
(since a^,j = ajUj > whenever rij > 0, see [5]), hence we have the partial fractions 
decomposition 

P n,N{x) ^(Z-Xfu/N' 
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where {x^i : 1 < i < \n\} are the zeros of and A^i > for every i < \n\. Let .fT be a 
compact set in C \ [0, oo), then ior x & K we have that 



Pn-g. N (x) 



\n\ 
i=l 



< 



X 



xn,i/N\ ~ 5 ^ 



^2 Ah >- 



where 



5 = inf { 1 2: -y\ ; z e K,y E [0, 00)} > 

is the minimal distance between K and [0, 00). Since P^n and Ph-s^n are monic poly 
nomials, one has X^i= 1 — 1, so that we have the bound 

Pn-e u N{%) 



Pn,N(x) 



1 



(3.2) 



uniformly for x G K. Take the recurrence relation (12. 4p with x replaced by Nx, and 
divide by C^(Nx), which is allowed since x G K cannot be a zero, then we find 



x 



Pfi+e k ,N( X ) , a k + \n 



+ 



PnAx) N 

If we use the bound (13. 2p . then this gives 



rijdj Pft-e^Nix) 



P n+ g ktN (x) _ a k + \n\ 



Ph,n( x ) 



N 



1 r 

< Ly 



3=1 



n j a j 

N 2 



Clearly, when n, N — > 00 in such a way that n/N — > t, we have 



lim 

n—>oo,n/N—>t N 



n n 

hm = £, 



and 



so that 



lim 



djTlj 



rij n 

....... , a,- — — - = 0, 

>oo,n/Jv^i iV 2 J n N 2 



lim 



n+e k ,N 



[x) 



uniformly for x G K, which proves the theorem. 



x — t, 



□ 



Observe that the same result will hold for any family of multiple orthogonal polyno- 
mials for which a^j > whenever rij > and 



n— >oo n 



lim — 



0. 



n— >oa n 



where rij = \_qjn\, with < qj < 1 and X)j=i = The ^ ac ^ that a n,j/ n2 ~^ simplifies 
the asymptotic analysis a lot and the limit function is an easy polynomial function of 
degree 1. In general, the asymptotic analysis for ratios of multiple orthogonal polynomials 
would involve a limit function which is the solution of an algebraic equation of degree r+1. 
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4 Asymptotic distribution of the zeros 

Next, we will obtain the asymptotic distribution of the (scaled) zeros of the multiple 
Charlier polynomials. For this, we introduce the zero measure 

1 \n\ 



i=l 



and we want to show that these (probability) measures converge weakly to a (probability) 
measure v t as n, N — > oo and n/N — > t > 0, which then describes the asymptotic 
distribution of the zeros. Again we will take multi-indices n such that rij = [nqj\ , where 
< q.j < 1 and Y7j=i 1j = 1> so that \n\/n — > 1 as n tends to infinity. In order to prove 
this weak convergence, we will investigate their Stieltjes transform 



dv n , N {y) 1 p k 



X 



x-y \n\ Ph,n( x ) 



iGC\[0, oo) 



where Pft^(x) = Cfi(Nx)/N^, and show that they converge to a function, which we can 
identify as the Stieltjes transform of a measure u t . The Grommer-Hamburger theorem [6] 
then tells us that the measures z/ n jy converge weakly to u t as n, N — > oo and n/N — > t. 

Theorem 4.1. Suppose rij = \_qjn\ , with < q,j < 1 and Y^j=i Qj = 1 an d ^at di > for 

1 < i < r and ai ^ aj whenever i ^ j . Let < x^,2 < ■ ■ ■ < xa,\n\ be the zeros of C^. 
Then 

1 N 1 /•* 

lim , mE fi^/N) = - / f{x) dx (4.1) 

n,N^oo,n/N^t \n\ z — ^ t J 

for every bounded continuous function on [0,oo). This means that the zeros of Crf(Nx) 
are asymptotically uniform on the interval [0,t] when n, N — > oo and n/N — > t > 0. 

Proof. We will prove that 

to 1^ = 1 /*_!_* (4.2) 

n,N^oo,n/N^t \n\ Pfi^{x) t J Q X — y 

uniformly for x G K, where K is a compact set in C \ [0, oo), which by the Grommer- 
Hamburger theorem (see, e.g., [6]) is equivalent with the weak convergence to the uniform 
measure on [0, t]. We will prove this by induction on r. For r = 1 we deal with the zeros of 
Charlier polynomials and the multi-index n is an integer which we denote by n. Observe 
that 

1 Pu,n( X ) 1 ( ' Pk+1,N\ X ) Pk,Ni x ) 



nP n , N { x ) n^\P k+liN (x) P k ,N(x) 
and straightforward calculus gives 

P Ul,N( X ) P Ln( X ) fP k +l,N( X )\' I fPk+lA X ) 



Pk+l,N{X) Pk,N{X) V P k,N{X) ) I V Pk,N{X) 

Hence we may write 

1 P u,n( x ) _ iyV Pk+i,N(x)\ j ( P k+1>N ( X 



nP ntN (x) n/ ^ \ Pk,N(x) J I \ P k ,N(x) 



We can rewrite the sum as an integral by putting k = \ns\ , so that 

1 Pti,n( X ) _ f l f Plns\ + l,N(x) \ j f P[ns}+l,N( x ) 



ds. 

nP n>N (x) J V Plns\,N(x) ) I \ P\ns\,JSl(x) 

Now we let n, N — > oo in such a way that n/N — y t, and we use Theorem 13.11 (with r = 1) 
to find that uniformly for x G K (K a compact set in C \ [0, oo)) 

i PnA x ) f 1 & - sty , 

lim — — — t-^t = / ds, 

n,N-*oa,n/N-+t 71 P n ,N\ x ) JO X — St 

where the ' in the integral is a derivative with respect to the variable x. The integral on 
the right is (use st = y) 

1 /"* 1 
ds — - I dy, 



x — st t J x — y 

which proves (14. 2p for r = 1. 

Now suppose that H4.2[) is true for r — 1. Observe that 

Pu,n( X ) Pn-n r e r ,N( X ) ST^ ( Pn-ke r ,N{ x ) \ / / Pn-ke r ,N( x ) \ , , .-, 



Pti,n( x ) Pn-n r e r ,N( x ) \ Pn-(k+l)e r ,N ( x ) J ' \Pn-(k+l)e r ,N( x ) 

The multiple orthogonal polynomial Prt-n r e r is in fact a multiple orthogonal polynomial 
with only r — 1 measures (/ii, . . . , /i r -i), hence we can use the induction hypothesis to find 

lim i fUfegW = i r {1 -"» _±_ iy (4 4) 

n,7V-s.oo |n| - n r P«_^.^,jv(ar) {1 - q r )t J x-y 

Note that (|n| — n r )/n —y 1 — q r , which explains the appearance of 1 — q r in the last 
formula. We can write the sum as an integral by taking k = \n r s\: 

1 ' Pn-ke r ,N( x ) V I ( P n-k)e r ,N( x ) 



r k=0 



Pn-(k+l)Sr,N\ x ) ) ' \Pn-(k+l)e r ,N( x ) , 

1 ( Pn-[n r s\e r ,N{. x ) \ if Pfi-[n r s\e r ,N{. x ) 



\Pn-(ln r s\+l)e r ,N(x) J ' \Pn-{\n r s\+l)e r ,N{x) 



ds. 



Now use Theorem 13.11 to find 

Pn-ke r ,N( x ) \ / / Pn-ke r ,N( x ) 



n r — l 

I 

lim 



h n r L — 4 

k=0 



n,N^oo,n/N^t 7l T \P n _( k+ x)g rjN (x) J I \Pn-(k+l)e r ,N( x ) 



' (x - (1 -*.).)• ^ = J_ _ iy 



x - (1 - g r s)t g r t ^(l-gr)t x ~ y 

where the last equality follows after using the substitution y = (1 — q r s)t. Note that 
(|n| — [n r s\)/n — y 1 — g r s, which explains the factor 1 — q r s in the asymptotic formula. 
Now combine (IP]) and (HT5)I in (jOjl to find 

1 PkM x ) l-Qr [ {1 ~ qr)t l,Or/" 1 . 
hm — — — — = — — / dy H / dy 



n,N^oo,n/N->t \n\ P H ,n( x ) (1 ~ Qr)t Jo X ~V Qrt J(l-q r )t X ~V 

— - I dy, 

t Jo x-y 

which proves (14. 2p . □ 
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5 Parameters depending on the degree 



We get more interesting asymptotics when some of the parameters depend on N and grow 
together with the degree \n\. The case where only one parameter depends on N can be 
worked out in detail. 

Theorem 5.1. Suppose rij = \_qjn\, with < qj < 1 and YTj=i1j = 1> so ^at \n\jn — > 1 
as n — >< oo. Consider Poisson distributions with parameters (ai, a2, . . . , a r -i, Na r ), i.e., 
the last parameter grows linearly with N. Then for t > one has 



lim 



C*n+e r (Nx) x — a r — t + — a r — t) 2 — Aa r q r t 



n,N->oo, n/N-H NCh(Nx) 

and for 1 < k < r 



g r (x) 



lim 



Cn+e k (Nx) 
n,N-»£>"n/N->t NCh(Nx) 

uniformly on compact sets of C \ [0, 00). 



x — t 



g r (x) 



(5.1) 



(5.2) 



Proof. We still use the notation Pn iN (x) = C^(Nx)/N n \ but now keep in mind that 
depends on the r parameters (a 1; . . . , a r _i, Na r ) so that the parameter N appears not only 
in the scaling of the variable (Nx) but also in the last parameter (Na r ). The recurrence 
relation (\2A\\ . after dividing by C^(Nx) gives for x G K, where K is a compact set in 
C\[0,oo), 



_ Pft +SkiN (x) a k + \n\ 



r-l 



Prt.Nlx) 



N 



+ nj<Xj - P "" g J- Af ( X ) + n r a r Pn-e r ,N(x) 



^ iV 2 P n>N (x) N P HtN (x) 



when 1 < k < r — 1, and for k = r we have 



_ ^n+ev,jv(z) iVa r + \n\ 



1 HjOj Pa-S h N{x) | n r a r Pfi-e r ,N(x) 



Ph,n( x ) 



N 



+ E 



- iV 2 P fJjAr (x) 



+ 



N P h ,n(x 



If we use (13.21). then for 1 < k < r — 1 



Pn+e k) N{x) 



Ph,n( x ) 



— X + 



a fc + |n| n r a r P^.^a^x) 



N 



+ 



r-l 



^E^ 



and for k = r 



so that 



P 



n+e r ,N 



X ) 



-Prj.vfa;) 



x + 



A^a r + Inl n r a r P^.^- ^(x) 



iV 



+ 



< 



1 r_1 
^E n i G i 



5iV 2 



i=i 



lim 

n—too,n/N—tt 



Pn +Sk , N (x) _ x + t+ t P*-erAx) 

PaNiX) 



and 



lim 

n— >oo,n/N— H 



Ph,n( x ) 

Pn+e r ,N(x) 



Kk<r-h 



Ph,n( x ) 



x + a r + t + a r q r t- 



Pn-e r ,N( x ) 



PnM x ) 
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uniformly on K. The bound (I3.2p implies that {P^g.^{x) / Ph,n(x) : n, N G N} is a 
normal family on every compact subset of C \ [0, oo), hence there is a subsequence which 
converges uniformly on K: 



and, by taking further subsequences, this convergence holds for every j for which 1 < j < 
r. With our previous estimates, this gives 

lim — i — J ' ' — x — t — a r q r th r (x), 1 < j < r — 1, (5.3) 

and 

hm — —^— = x — a T — t — a r q r th T {x). (5.4) 



A technical estimation (see Lemma 15.11 at the end of this section) implies that 

Pfi,N{x) PfL-e^Nix) 



lim 

n— too, n/N—tt 



Pn+ej,N\ x ) Pn,N\X 

uniformly on K, hence (\5A\\ gives 

1 



0, 1 < j < r, 



h r (x) 

If we put g T {x) = l/h r (x), then this gives a quadratic equation for g r (x), with solutions 

x — a r — t ± \J(x — a r — t) 2 — 4a r q r t 
2 ' 

Since h r (x) = 1/x + 0(l/x 2 ), we need to choose the solution with the positive sign 
for g r {x). This limit is independent of the subsequence that we selected, hence every 
convergent subsequence has the same limit, which implies that the full sequence converges 
to this limit. This gives (I5.ip . and by using (15. 3 p we easily find (15. 2p . 

□ 

The limit function g r (x) is the solution of a quadratic equation. In general, if k < r 
of the parameters grow linearly with N, then the limit function is expected to be the 
solution, which grows as x when x — > oo, of an algebraic equation of degree k + 1. 

For the asymptotic behavior of the zeros we have 

Theorem 5.2. Suppose rij = \_qjn\, with < qj < 1 and Y^j=i Qj = ^ so that \n\/n — > 1 
if n — > oo. Consider Poisson distributions with parameters (a±, a%, . . . , a r _i, Na r ), i.e., 
the last parameter grows linearly with N . Then for t > one has 

I W i r(i-9r)t rPt 

lim T^/J{xti,j/N) = - f(x)dx + q r v(x)f(x)dx 

n,N-^oo, n/N-+t \7l\ ^ t J Q J at 

for every bounded continuous function f on [0, oo), where (1 — q r )t < a t < (3 t and v is a 
probability density on [a t ,Pt\- 
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Proof. We can start from equation (14. 3ft : 



Pn,N\X) Pfi^ nr g ri N(x) \Pn~(k+l)e r ,N(x) J I \ Pfi-(k+l)e r ,N (x) 

The multiple orthogonal polynomial Pn- nr g r is in fact the multiple Charlier polynomial 
with the r — 1 parameters (ai, . . . , a r -i), which do not depend on N. Hence we can use 
Theorem 14.11 which gives (14. 4p . We write the sum as an integral, as we did in the proof 
of Theorem 14.11 but now we use Theorem 15.11 to find 

hm ly ( P ^n(x) V / ( Pn-^N{x) \ = fgrjx,*)' 

n,N^oc,n/N^t U r ^ \ Pfi-( k+1 )g rtN (x) J I \ Pn-(k+l)e r ,N {x) J J 9r{X, S) 

where 

x — a r — (1 — q r s)t + y/Tx — a r — (1 — q r s)t) 2 — Aa r q,Al — s)t 
9r(x, s) = 

and the prime is the derivative d/dx. The g r (x, s) is obtained from Theorem 15. II after the 
substitutions 

Qj ->• : 1 < J < r - 1, q r -> ^ — , 

1 — g r s 1 — g r s 

so that £^ =1 <?j = 1, 

n r — > n r — \n r s\ , n — > n(l — q r s), t — > (1 — q r s)n. 

Observe that 

9r(x,s)' 1 



fivO, s) - a r - (1 - g r s)t) 2 - Aa r q r (l - sji 

and if we use the well known Stieltjes transform 

1 If 1 ! dy 



[ — 7=, *eC\[-i,l], 

J-ix-y J\ - v 2 



Vx 2 - 1 TT J^X-y y/l - yl 

then one finds 

g r (x, S y 1 1 



c/ r (x,s) it J a(s) x-y ^Aa r q r {l - s)t - (y - a r - (1 - q r s)t) 2 ' 

where 

a(s) = a r + (1 — g r s)t — 2 a/ a T .g r (l — s)£, /3(s) = a r + (1 — q r s)t + 2\/a r g r (l — s)t. 
In order to write 

f l 9_r{x i Sy ds 
Jo 9r{X,S) 

as a Stieltjes transform, we need to change the order of integration in 
- 1 1 f^ s) 1 dy 



o * Ja(s) x ~y xjAa r q r {\ - s)t - (y - a r . - (1 - q r s)t)' 

11 



ds. (5.5) 



Observe that 

a(0) = a r + t — 2y / a r q r t, (3(0) = a r + t + 2^a r q r t, a(i) = (3(1) = a r + (1 — q r )t, 

and that the function (3 is monotonically decreasing for s G [0, 1]. We need to distinguish 
between two cases. 

Case 1: a r > q r t. In this case the function a is monotonically increasing for s 6 [0, 1], 
see Figured! 

3 

P 

2.5- 

2- ^\ 



0.2 



0.6 



Figure 1: The functions a and (3 for case 1 



If we define a t = a(0) and (3 t = (3(0) then 



-1/ \ -y-a r +t+2yj a r (y-(l-q r )t) , , / n n, 

7-1, ■ -y-a r +t+2y/a r (y-(l-qr)t) . /-, \ , ^ ^ n 

P (y) = ^ , «r + (i - Q r )t <y < Pt 

so that interchanging the order of integration in (I5.5P gives 

ds 



-y — ar + t + 2^/ a r (y— (1 -qrjt) 



1 F Pt dy 

7T Ja t x-yJo y/4a r q r (l - s)t - (y - a r - (1 - q r s)t) 2 

When we change the variable s to a new variable u by 



-y - a r + t + 2uy/a T (y - (1 - q r )t) 



q r t 



then the integral simplies to 



nq r t 



dy 



du 



2yJa r (y-{\-q r )t) 
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This gives the weight function 



v{y) 



l 



irq r t 



du 



y + a-r — t 



a t <y< 0t- 



2^a r (y-(l-q r )t) 

An easy exercise gives that (1 — q r )t < a t < (3 t . 



Case 2: a r < q r t. In this case a has a global minimum on ]0, 1[ at s = 1 — a r /q r t, and 
the minimum is (1 — q r )t, see Figure |2j 



2.8- 



2.6- 



2.4- 



2.2- 



2- 



1.6 



1.4- 




0.2 



0.4 



0.6 



Figure 2: The functions a and (3 for case 2 



Interchanging the order of the integrals in (I5.5P now gives two pieces 

-V — a r + t + 2y' a r (y — (l — q r )t) 



+ 



1 r (0) _dy_ 

K J(l-q r )t X - y J -y- a r+t-2Var(y-(l- q r)t) ^O^^l - s)t - (?/ - <2 r ~ (1 - q r S)t) 2 

1 fPV) dy 



qrt 



qrt 

-y — ar + t + 2yj ar(y — (l — qr)t) 

qrt 



ds 



ds 



* Ja(o) x ~yJo ^/Aa r q r (l - s)t - (y - a r - (1 - q r s)t) 2 ' 

The change of variable s — > u with 



— y — a r + t + 2u^a r (y — (1 — q r )t) 
q r t 



now gives 



I r (0) dy f 1 du 1 f m dy f 1 



*q r t J(i- qr )t x-y J-i y/l - u 2 irq r t J a(0) x-y J v+* r -t \/Y^ 



du 



2^/a r (y-(l-qr)t) 



So if we now define a t — (1 — q r )t and (3 t = (3(0), then obviously (1 — q r )t = a t < (3 t 
and the weight function becomes 

v(y) = — + , (1 - q r )t <y< a(0) 
q r t 
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and 

v(y) = ^~f ~J^, a(0) < y < 0(0). 

2y/a r {y-{\-q r )t) 

So in both cases we get 

n,N^oo,n/N^tn r \P H _ (k+1) g r:N (x)J I \ P H _ ( k +l)e r ,N (x) J J at X-y 

and combining this with (14.4D gives 

r 1 1 [ {1 - Qr)t 1 . ^ [* v(y) 
hm — — ■ — — = - / ay + q r / ay, 

n,N-^oc,n/N->t \n\ Pn,N{x) Uo X — y J a t X ~ V 

which gives the desired result in view of the Grommer-Hamburger theorem [6]. □ 

The first portion of (1 — q r )n of the zeros of C^{Nx) are uniformly distributed on 
[0, (1 — q r )t] and hence the constraint that 'between two positive integers there can be 
at most one zero' is in action and the zeros are forced to approach the first (1 — q r )n 
integers in N. If a r > q r t (case 1) then the last portion of q r n of the zeros have a different 
distribution on an interval [a t , f3 t ] = [a r + t — 2y/a r q r t, a r + t + 2y/a r q r t\ to the right of 
the interval [0, (1 — q r )t] where the other zeros accumulate. This means that those last 
q r n zeros are less dense distributed and some of the intervals between two integers may 
be free of zeros. If a r < q r t (case 2) then some of the q r n last zeros are still uniformly 
distributed on [(1— q r )t, a r +t — 2y/a r q r t\ but the remaining zeros are less dense distributed 
on [a r + 1 — 2y/a r q r t, a r + 1 + 2y/a r q r t\ and this interval now touches the interval where the 
zeros are uniformly distributed. In fact, a transition occurs when a r = q r t in the sense 
that the (scaled) zeros have a zero distribution on two disjoint intervals when a r > q r t 
and the zero distribution is supported on one interval when a r < q r t. Moreover, since 

1 du 



CVl - z, z->l- 



u- 



and for a r > q r t 



y + a r -t JCi(y-a(0)), y -)> a(0)+ 



2y/ar(y-(l-q r )t) \C 2 (f3(Q) - y), y ^ 0(0)- 

we see that the density v near the endpoints a(0) and /3(0) tends to zero as sjy — a(0) 
and a//3(0) — y, respectively (see Figure [3], picture on the left, for a r = 1, q r = 1/10 and 

For a r < q r t we still have v(y) ~ Cy/P(Q) — y near the endpoint 0(0). The transition 
from uniform density to non- uniform density occurs at y — a(0), but now 

y + a r — t 

->■ -1, y^a(0) + 



2y/a r (y - (1 - q r )t) 



so that v(y) — > l/q r t as y — > a(0)+, and the density is continuous at the transition point 
a(0) (see FigureEJ picture on the right, for a r = 1/10, q r = 1/5 and t = 1). 
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Figure 3: The density of the (scaled) zeros of multiple Charlier polynomials 



When a r = q r t we have 

y + a r — t 



2^/a r {y - (1 - q r )t) 



->0, y^a(0) + 



so that v(y) — > l/2q r t as y — > a(0)+, so that the density is not continuous at the transition 
point. 

Such transitions also occur when k < r of the parameters depend linearly on N. In 
that case the zeros of Ca(Nx) may accumulate on at most k + 1 disjoint intervals. If all 
the parameters depend on N (i.e., k = r) then the zeros accumulate on at most r disjoint 
intervals. The analysis for k > 1 is more involved since this involves algebraic functions 
of order k + 1. 

One technical, but crucial, step in the proof of Theorem 15.11 is the following. 

Lemma 5.1. Let P^ N {x) = Cn(Nx)/N^, whereC^ are the multiple Charlier polynomials 
with parameters ( Na r ). Let K be a compact set in C\ [0, oo) ; then for every 

k and £ with 1 < k, t < r one has, uniformly for x G K 



lim 



0. 



Pn+e k ,N(%) Pn+e k -e e ,N(%) 

Proof. From the recurrence relation, we have 

_ Pfj+g kiN (x) a k + \n\ ^ njcij Pfi-g hN {x) n r a r Ph-s^n{x) 

X ~ Ph,n(x) N jr[ N 2 P n>N (x) N P HtN (x) 

when 1 < k < r — 1, and for k = r we have 

_ Pfj+g r ,N{x) Na r + |nj njaj Pn-g^Njx) n r a r Pn-g r ,N{x) 

X ~ Pu,n{x) N N* P n , N (x) N P^ N (x) 

We will denote 

Pn-g^Njx) 
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and the bound (13. 2p then gives 



r-l 



\En(x)\ < 



1 C\n 



< 



3=1 



SN 2 



where C > is a constant (in fact on may take maxi<j< r _i aj). If we change n to n — 
then 



Pft + g k -g hN (x) a k + \n\ - 1 



Pn-ee,N\X) 



+ 



+ E H _ s .(x) 



,n r — 5 r i)a r Pfi_ 



■e r —et 



N 



Pn-e f ,N\X) 



when 1 < k < r — 1, and for k = r 



P 



x 



n+e r —§£,N 



'x) Na r + \n\ — 1 

+ - 



Efi—gf (x^j -\ 



(n r — 8 r p)a r Pi 



r u r,l) KJj r ± n—e r —ei 



,N[X) 



Pn-etAx) ' N ' " Ce " y ' N Pn-e t Ax) 

If we subtract the equations for n — ee from those with n, then we find 

Pn+e k ,N( x ) Pri+e k -e e ,N{ x ) 1 







PnM X ) Pn-g f A X ) N 



+ 



n r a r i Pf[—g r} N{x) Pfi-e r -e t ,N(x)\ 5 r /a r Pfi-g r -g e ^]\r(x) 



N V Pn, N (x 



We have 



\Eft(x) - Eft-g f (x)\ < 



Pfi—et,N\Xj 

2C\n\ 



+ 



N P 



n—ee,N 



x) 



and we will take \n\ < C 2 N, therefore we have 



Pn+e k ,N(x) Pfi + g k -g e ^(x) 



Ph,n(x) 



Pn-en,N\X) 



< ^ + C 2 a r 



6N 2 



Pn-e r ,N^X) Pn-e r -e e ,N(x) 



Pn,N(X) 



Pn-e e ,N{x) 



where C\ and C 2 are constants. If we use the bound (13. 2L then 



Pn+e k ,N(x) Pri+g k -e t ,N\X) 



P*..n{x) 



P 



n—ee,N 



A 



>s 2 



Pn,N(x) Pfi-e e ,N(x) 



so that 



Ph,n{x) 



Pn-e p ,N(x) 



\Pn+e k ,N(x) Pf£ + g k _g e ^(x) 

If we use the notation 



^ C\ C 2 a r 



Pft+e k ,N(x) Pfi^g k -g l ^(^x) 

Pfi-e r ,N(x) Pfl-g r -g l ^( y X) 



N5 3 5 2 



D 



Prt.N(x) 



Ph,n{x) 

Pn-ep,N{x) 



Pn-e f ,N\X) 



Pn+e k ,N(x) Pti + g k _g tj ]y(x) 



then this gives 



Put 



Ci C 2 a r n 



Dm = max D H k e , 

Kk<r 
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then one has 



Am < 



+ 



C 2 a, 




N5 3 



Iterating this inequality gives 




3 



Now choose a compact K' (with an accumulation point) far enough from [0, oo) so that S 
is large and C^Orjb 2 < 1. Then for x G K' 



The bound (13. 2p gives < 2/5, hence if we put n = ( |_^<?ij , • • • , [nq r \ ) and let n,N 
oo such that n/N — > t > 0, then 



uniformly for x £ K' . So we have convergence of D^^ — > uniformly on a set K' with 
an accumulation point, but then Vitali's theorem implies that converges to zero 
uniformly on every compact K where a bound (13 .2\\ holds, hence for K C C\ [0, oo). □ 

6 Concluding remarks 

In this paper we have investigated the ratio asymptotic behavior of the multiple Char- 
lier polynomials and from it we obtained the asymptotic distribution of the zeros, after 
proper rescaling. The next step is to find the asymptotic behavior of the polynomials 
themselves: the strong asymptotic behavior or the uniform asymptotic behavior. As in 
the case of the usual Charlier polynomials, one will need to look at different regions in 
the complex plane: away from the positive real line, on the oscillatory region where all 
the zeros are, near the largest zero, near the origin, etc. One way to do this is to use the 
integral relation which can be obtained from the multivariate generating function and to 
apply a steepest descent analysis (but for a multiple integral), as was done by Goh [7] and 
Rui and Wong [20] for Charlier polynomials. Another way is to use the Riemann-Hilbert 
problem (for (r + 1) x (r + 1) matrices) and the steepest descent method for oscillatory 
Riemann-Hilbert problems, as was done by Ou and Wong [17] for Charlier polynomials. 
One of the steps in that asymptotic analysis is to transform the Riemann-Hilbert prob- 
lem to a normalized (at infinity) Riemann-Hilbert problem, and this requires ^-functions 
which are logarithmic potentials of the asymptotic zero distribution. Hence the results in 
Section 4 and Section 5 (in particular Theorem 15. 2j) will be needed. 
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